Reducing nonideal to ideal coupling in random matrix description of 
chaotic scattering: Application to the time-delay problem 
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We write explicitly a transformation of the scattering 
phases reducing the problem of quantum chaotic scattering 
for systems with M statistically equivalent channels at non- 
ideal coupling to that for ideal coupling. Unfolding the phases 
by their local density leads to universality of their local fluctu- 
ations for large M. A relation between the partial time delays 
and diagonal matrix elements of the Wigner-Smith matrix is 
revealed for ideal coupling. This helped us in deriving the 
joint probability distribution of partial time delays and the 
distribution of the Wigner time delay. 

PACS numbers: 05.45.-a, 24.60.-k, 73.23.-b 



The random matrix theory (RMT) is generally ac- 
cepted to be an adequate tool for describing various 
universal statistical properties of quantum systems with 
chaotic intrinsic dynamics, see Ref. and references 
therein. In particular, one can distinguish two variants 
of the RMT approach allowing one to address the chaotic 
nature of quantum scattering. The first one considers 
the scattering matrix S as the prime object without any 
reference to the system Hamiltonian. The probability 
distribution P{S) of S at the fixed energy E of incident 
particles is chosen to satisfy a maximum entropy principle 
and natural constraints which follow from the unitarity 
and causality of S, and the presence (or absence) of the 
time-reversal (TRS) and spin-rotation (SRS) symmetries, 
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Such a distribution is known as the Poisson kernel 
and uses the phenomenological average (or optical) S'- 
matrix S(E) as the set of input parameters. Without 
loss of generality S can be considered as diagonal j|. 
P(S) depends also on the number of scattering channels 
M and the symmetry index (3 [P=2 for a system with 
broken TRS, and /3=1(4) if the TRS is preserved and the 
SRS is present (absent)]. 

The approach proved to be a success for extracting 
many characteristics important in the theory of meso- 
scopic transport However, correlation properties of 
the S'-matrix at close values of energy E as well as spec- 
tral characteristics of an open system related to the so- 
called resonances turn out to be inaccessible in the frame- 
work of such an approach, essentially because of the one- 



energy nature of the latter. To address such quantities 
one needs to consider the Hamiltonian H of the quantum 
chaotic system as the prime building block of the the- 
ory. It amounts to treating if as a large N x N random 
matrix of appropriate symmetry and relating S to the 
Hamiltonian by means of standard tools of the scatter- 
ing theory ||^] . This idea supplemented with the super- 
symmetry technique of ensemble averaging || resulted 
in advance in calculating S-matrix correlation functions 
|^|,^| and many other related characteristics such as, e.g., 
time delays Jlp|-|T^|, see Refs. |pr|,pj for a review. 

In the limit N—>oo one can prove |Ts)| the equivalence 
of both mentioned approaches by deriving the Poisson 
kernel ([!]) from the Hamiltonian approach (see also Ref. 
|]llf ) , with the average S-matrix being 



S(E) ab 



l--y a [iE/2 + w{E)] 
l + f a [iE/2 + irv{E)] 
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Here, the average density of states 
(E/2) 2 determines the mean level spac- 



indcpendent of (3 
v(E) = n- 1 ^/T~ 

ing A=(^A^) _1 of the closed system, and phenomenolog- 
ical constants 7 C > characterize the coupling strength to 
continuum in different scattering channels (c=l,. . . ,M). 

The particular case of ideal coupling, S=0 [when the 
transmission coefficients equal unity for all channels, see 
Eq. (||) below] , plays an especiallly important role for the 
S-matrix approach |^| . Equation ([!]) simplifies then to 
Po(S)=const, which is invariant under the transforma- 
tions of S, leaving the measure invariant. Such a situ- 
ation corresponds to the so-called Dyson's circular en- 
semble (CE) of unitary matrices and is much simpler to 
handle analytically. 

A general situation of nonideal coupling, S^0, turns 
out to be much more complicated. It is natural to ex- 
pect, however, that results obtained for the case of non- 
ideal coupling could be related to those at ideal coupling. 
Although many useful ideas around such a relation were 
discussed in the literature [^|,[l3 11 ll| , we are not aware 
of explicit relations, to the best of our knowledge. 

In this Rapid Communication we consider the most 
simple but physically important case of statistically 
equivalent channels. We demonstrate the validity of 
the following simple statement (and discuss several ap- 
plications of it): Let S(E) = Us(E)U\ where s(E) 



diag(e 



2iS 1 (E) 



>e 2lS M(E)^ kg ran( j om S'-matrix at 



the energy E, the distribution of which is given by the 



1 



Poisson kernel dl|) with an explicit parameterization of 
S(E) from Eq. (g) (7 C =7 for all c). Then for every E the 
transformation of the eigenphases S c (E) 



arctan 



1/1 E 
— 7—7 ( -tanSJE) + — 
ttv{E) V7 2 



(3) 



maps them to the eigenphases <p c of the random scatter- 
ing matrix, the distribution of which is given by the CE of 
the same symmetry. In particular, the joint probability 
density function (JPDF) of cf> c is § 



Po({0c}) oc Yl 



3 2#„ 



= 2# b I 



(4) 



a<b 



The matrix U of energy dependent eigenvectors uni- 
formly distributed in the orthogonal, unitary, or symplec- 
tic group (for (3=1,2, or 4, respectively) is not affected by 
map (H) . This is a consequence of statistical equivalence 
of the scattering channels. 

The suggested transformation was first noticed and 
verified in Ref. jn) for the case of broken TRS ((3=2) 
and further exploited in Ref. jl2| . It can be easily gener- 
alized to the other symmetry classes as follows. We cal- 
culate first the Jacobian of the transformation (||) . After 
a simple algebra it can be represented as 



cos" 



T 
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where 



T(E) = 1 - \S{E)Y = 



4jnu(E) 



1 +7 2 + 2-/ttis(E) 
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(6) 



is the energy dependent transmission coefficient Q. We 
note that exactly the same factors as Eq. (||) appear in 
Eq. (|l|) . Employing the identity 



e 2i0 6 | _ \2iS a 



COS (p a COS <fib 



irvj cos 8 a cos 5b 
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we substitute Jacobians (||) into Eq. (||) and, making use 



of the identity ]jf b f a f b = (J] c fc) M \ arrive at 



P({S C }) oc [] 
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With Eq. (|^) taken into account, we immediately rec- 
ognize in this expression the JPDF of the eigenphases 
corresponding to Poisson kernel ([!]). Due to the scalar 
nature of transformation (Q) it docs not change the ma- 
trix U of eigenvectors. 

Let us start with considering the mean density p(5) of 
scattering (eigen)phases at arbitrary coupling. It is self- 
evident that the phases in the CE (i.e. for the case S=0) 
are uniformly distributed on the unit circle, the average 
density being merely p (</>) = 0-/M) ^(4>-4>c) = 1/tt. 



The corresponding density for Sy^O is not constant. In- 
deed, using the identity p(S)dd = po(4>)d((i, we see that 



p(S) = - 

7T 



dj)_ 
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Although simple, this relation is an important one and es- 
tablishes the physical meaning of the Jacobians of trans- 
formation (|J) relating them to the corresponding densi- 
ties of the scattering phases. Density ([)]), being expressed 
in terms of S only, does not depend on the particular 
choice of S used in the derivation as long as the average 
S'-matrix is proportional to the unit matrix. 

It is instructive to look at Eq. (||) in the limit of large 
number of channels when the typical difference S a — 6b ~ 
1/M <C 1. Then one can expand 5 C = So + S c (<5 c <Cl) 
around, say, 5q. The leading contribution is given by 

P(R})°C Ua<Ja-J b \ Uc\^o/d5/^ M - 1)/2+ \ Which 

further goes to Po({0c})ocn a <6 \4>a 4>b\^ and agrees with 
distribution ([|) of the CE upon the proper rescaling of 
the phases, 



4>c = \d(pc/dS c \ s S c = irp(S )5c 



(10) 



We see that in the limit M 3> 1 the local fluctuations 
of the phases unfolded by their local density turn out to 
be uniformly described by the CE at arbitrary coupling 
strength. Such a universality in statistics of phases of 
random unitary (scattering) matrices has much in com- 
mon with that typical for eigenvalues of random Hamil- 
tonian matrices []l| and is in agreement with results of 
realistic numerical simulations for M = 23 [ flGf . 

Let us now consider an application of the same ideas to 
the time-delay problem, where such a universality reveals 
itself explicitly. Following the original wave-packet anal- 
ysis by Eisenbud, Wigner and Smith Jl?]] it is natural to 
define [|lTj| the partial time delays via the energy deriva- 
tive of the scattering phases, t c = 2TidS c /dE. Their sta- 
tistical properties have been studied in much detail in 
the framework of the Hamiltonian approach for the case 
of broken |ll[] and preserved TRS as well as in the whole 
crossover region of gradually broken TRS [[l2| . Recently, 
some of these predictions were successfully verified on the 
model of a quantum Bloch particle chaotically moving in 
a superposition of ac and dc fields JT^ . 

In particular, the mean density of partial time delays 
7 , (r) = (1/M) J2 C ^i T ~ T c) turns out to be especially sim- 
ple at ideal coupling, T=l, when it reads as 



V (t = T/t H ) 



((3/2f M l 2 e-P' 21 

T((3M/2) i/3M/2+2 



(11) 



with tH=2TTh/ A being the Heisenberg time. Due to 
Eq. (|Io|), the partial time delays at ideal and nonideal 
coupling (tc°^ and r c , respectively) are simply related as 

ri°) = 2Wct> c /dE = np(S c )Tc . (12) 
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Here, we have neglected the smooth nonresonant depen- 
dence of p(5) on E. Since the phase and its derivative 
(the partial time delay) are uncorrelated quantities in the 
CE [^9|, their joint distribution factorizes: po{<j), t^) = 
{1/tv)Vq{t^). This is not the case for p(S,r), when 
Sy 0. The relation p (4>,T^)d^dT^ =p(S,r)d5 dr be- 
tween them allows us, however, to represent the density 
of partial time delays at nonideal coupling as 



a(#(j),T<°>(j,T)) 



Jo 9(5, t) 
= r^[np(S)} 2 V (7Tp(6)r). (13) 

Jo 7T 

One can easily convince oneself that such a formula 
reproduces in every detail the expression obtained in Ref. 
Jll[ by means of supersymmetry calculations. It is worth 
mentioning that the density of phases (|^) is independent 
of the underlying symmetry and therefore Eq. ( p"3| ) is also 
valid for the crossover regime of partly broken TRS. [Note 
that in the crossover regime Vo(t) is a slightly more com- 
plicated function, see Ref. 

Expression ( |l3| ) is the proper one for generalization to 
the JPDF of the partial time delays, w({t c }). Before 
doing this, we first establish a useful relation between 
t c and the matrix elements of the Wigner-Smith time- 
delay matrix Q = ~ih{dS/dE)S^ @- Writing S in the 
eigenbasis representation as S = UsU\ one obtains 



U^QU 



-%h 



ds 
"BE 1 



ih 



dU 
~dE 



(14) 



where [, ] denotes the commutator. The matrix s being 
diagonal, the diagonal elements of the second term in 
Eq. ( |l4| ) are zero, whereas the first term is exactly the 
diagonal matrix of the partial time delays. Thus, the 
partial time delays coincide with the diagonal elements 
of the time-delay matrix taken in the eigenbasis of the 
scattering matrix, 



r c = [U^QU] 



(15) 



The physical meaning of the diagonal elements of the 
time-delay matrix is well known: they describe the time 
delay of a wave-packet incident in a given channel @|lj- 
Thus, relation ([ll]) sheds more light on the physical 
meaning of the somewhat formally defined partial time 
delays. In particular, one expects that for the case of 
ideal coupling the inherent rotational invariance of the 
problem makes all the basises statistically equivalent and 
thus the JPDF of diagonal elements of the Q-matrix 
should coincide with that of partial time delays. 

The latter claim can be substantiated as follows. Fol- 
lowing the insightful paper [|9j, it is convenient to con- 
sider the "symmetrized" time-delay matrix Q s 



-ihS-W^S- 1 ' 3 . (16) 



This similarity transformation unveils the symmetry 
which is hidden in Q: Q s is already a real symmetric 
(hermitian, or quaternion self-dual) matrix for (3 = 1, (2, 
or 4). In the eigenbasis of S the diagonal elements of Q s 
and those of Q coincide. Moreover, in the case of chaotic 
scattering with ideal coupling, the matrix Q s turns out to 
be statistically independent of S, their joint probability 
density being P (S,Q S ) = P (S)W {Q S ), where 



W (Q S ) oc fl(Q s ) det(Q s ) 



-3 / 3M/2-2+^ e -(/3/2)t ff trC 



(17) 



is the probability density of the time-delay matrix p9[ . 
The latter is manifestly invariant under the choice of the 
basis for Q s proving the above statement on the relation 
between statistics of partial time delays and diagonal el- 
ements of the Wigner-Smith matrix. 

To find the corresponding JPDF u>o({r c }) one has to 
integrate out all off-diagonal elements of Q s which is a 
hard problem in general. For the case of unitary sym- 
metry, = 2, one can perform the job by splitting the 
integration into that over the matrix q = diag(gi, q^) 
of eigenvalues of Q s and that of the eigenvectors, V, 



<({Tc}) 



det (q) 



3M 



2({r c }), (18) 



with A(q) — Y\ a<b {q a — Qb) being the Vandermonde de- 
terminant. Here Q({r c }) = Jd[V] Uc=i 6 ( T c- (VqV j )cc) 
stands for the remaining integral over the unitary group 
which can be done, following Re f. p2[ ], by means of the 
famous Itzykson-Zuber formula |23[ . Finally, we find it 
more convenient to define the generating function of par- 
tial time delays rather than the JPDF itself and obtain 



-i(k 1 T 1 + ...+k M T M ) 



det [V>j(fc0] 

Ua<b( k a - fa) ' 



(19) 



where ipjQei) = f™dq q j-3M e -ik,q-t H /q^ the index j 
spans the values I = 1, . . . , M and j = 0, 1, . . . , M — 1. 

Such an expression allows us to calculate all the mo- 
ments and correlation functions of partial time delays by 
a simple differentiation. Moreover, setting in the preced- 
ing equation k\ = . . . =kM=k and calculating the corre- 
sponding limit in the right-hand side, we come to a con- 
venient representation for the distribution 'P]^(t w ) of the 
Wigner time delay, t w = (ri+...+TM)/MtH, for a system 
with broken TRS and ideal coupling to continuum, 



where ■tp < - n) (k) = d n ip :j (k)/dk n , and j, n = 0, . . . , M-l. 

The distribution of the Wigner time delay was earlier 
calculated explicitly only for the case of M=l |Tl] , |T2| ,24 , 
when it follows from Eq. (pi]). Compact expression (2C) 
is valid for /3=2 and arbitrary M j2f|. For M=2, Eq. (2C) 
can be integrated further to yield 



oc 



dke 



iMkt u 



det[^ n) (fc)] 



(20) 



3 



V 2 (t w ) oc t-W+Ve-^Um±, 2/3 + 2, £) , (21) 

with [/(a, 6, z) = [l/r(a)]/ ° o dyy a - 1 (l + y) b - a - 1 e~ 2y be- 
ing the confluent hypergeometrical function. Here we 
represented the above distribution ( pl| ) in a form covering 
all f3 — 1,2,4 which will be verified below. In particular, 
the asymptotic behavior at t w » 1 is V%(t w ) oc t~^~ 2 
in agreement with the known universal tail f-P M / 2 - 2 ( 
which is typical for the time-delay distributions in open 
chaotic systems p|,^2| p|JT^] . 

To verify Eq. (|21|) for /3=1,4, it is convenient to con- 
sider a general problem of finding the distribution W (Q) 
of the n x n submatrix Q standing on the main diagonal 
of Q s . This distribution is found to be 



VF (Q)oc#(Q)det(Q) 



-P(M/2+n-l)-2 -ptHteQ /2 



(22) 



The particular case n=l reproduces result ( |Tl| ) of the 
Hamiltonian approach. Equation ( [22] ) for n=2 helps in 
calculating the joint distribution Wo(ti,t2) of two partial 
time delays t\$ — Ti^/tjj for arbitrary M. One obtains 



w Q (t 1 ,t 2 ) 
Po(*i)^b(*a 



(23) 



The knowledge of ^(ii, £2) allows us to find further the 
distribution of the Wigner time delay for M=2 and thus 
prove the formula ([H]) for any (3. As follows also from 
Eq. (p3f) , there exist nonvanishing correlations between 
the partial time delays. They are, however, of differ- 
ent nature as compared to the correlations between the 
proper time delays (the eigenvalues of Q) which show re- 
pulsion p9| . 

For 5^0 the matrices S and Q s cease to be statis- 
tically independent variables and do correlate. There- 
fore statistical properties of diagonal elements of Q in 
arbitrary basis (save the eigenbasis of S) are different 
from that of partial time delays, unless coupling is ideal. 
Still, the JPDF u>({r c }) of the partial time delays at 
nonidcal coupling can be found by repeating basically 
the same steps which lead to Eq. (|l3|). The identity 
p({S c },{r c })d[S}d[r} = p o ({0 c },{ri O) })dMd[r(°)], to- 

gether with the statistical independence of </> c and tJ ^ 
[which follows from Eq. ([l7j)], allows us to relate w({t c }) 
and iuo({r c }) as follows: 



W ({r c }) = d[S] 



3({<M,{r c (0) }) 



d{{S c },{r c }) 



Pb({<M>te (0) }) 



M 



= dlS}Hlnp(5 c )}p({S c })w ({np(6 c )T c }), (24) 



where d[8] means the product of differentials. 

In conclusion, we suggest the transformation of the 
scattering phases, allowing one to reduce the problem of 



quantum chaotic scattering with statistically equivalent 
channels at arbitrary coupling to that for ideal coupling. 
Applications of this transformation to statistical proper- 
ties of phases and those of time delays are discussed. 
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